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We present an analysis of the mass of the X(3872) reconstructed via its decay to J/iIjtx^-k~ using 
2.4 fb~^ of integrated luminosity from pp collisions at -^i = 1.96 TeV, collected with the CDF II 
detector at the Fermilab Tevatron. The possible existence of two nearby mass states is investigated. 
Within the limits of our experimental resolution the data are consistent with a single state, and 
having no evidence for two states we set upper limits on the mass difference between two hypothetical 
states for different assumed ratios of contributions to the observed peak. For equal contributions, 
the 95% confidence level upper limit on the mass difference is 3.6 MeV/c^. Under the single-state 
model the X(3872) mass is measured to be 3871.61 ± 0.16 (stat) ± 0.19 (syst) MeV/c^, which is the 
most precise determination to date. 
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The discovery of the X(3872) [H, Q and many addi- 
tional unexpected states [3| has revived general interest 
in spectroscopy in the charmonium mass region. Initial 
attempts to explain the X{3872) as a conventional bomrd 
state of a c-quark and an anti-c-quark have shortcomings 
Q which triggered the development of unconventional 
explanations. Two popular models are a molecular state 
composed of 1?° and mesons [1, 0] , and a four-quark 
state 0. 

In efforts to resolve the nature of the X{3872), several 
of its properties have been measured. The first deter- 
minations of its mass [H, 0, S, 0| resulted in values very 
close to the D'^D*^ mass threshold. The observed width 
in these measurements was compatible with zero. Stud- 
ies of the X(3872) production properties in pp collisions 
suggest that the production mechanisms are sim- 
ilar to those for the '0(25') charmonium state. Several 
measurements constrained the quantum numbers spin 
(J), parity (P) and charge- conjugation parity (C) of the 
X(3872). These include evidence for the decay modes 
X(3872) ^ J/V'7, J/ipuj, and ^{2S)-f [lH, and a mea- 
surement of the mass distribution of the dipions from 
the X(3872) — > J/iPtt~^tt~ decay These measure- 

ments indicate an even C parity. A subsequent angular 
analysis constrained the quantum numbers to only two 
possibilities, J^^ = or 2 A possible further 

decay mode of the A'(3872) was identified as a peak near 
threshold in the D^D'^tt^ invariant mass spectrum 
with a mean mass more than 3 MeV/c^ above measure- 
ments in the J/07r+7r~ mode. Despite efforts on both 
the experimental and theoretical sides, the nature of the 
X(3872) still remains an unresolved puzzle. 

A measurement of the AT (3872) mass with increased 
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precision can provide crucial information for understand- 
ing its nature. Under the hypothesis of a molecular state 
the mass of the X(3872) has to be lower than the sum 
of the and D*^ masses. The four-quark state hy- 
pothesis predicts the existence of two distinct particles 
that differ by the light-quark content bound to the cc 
quarks. These two particles should have slightly differ- 
ent masses, and the model of Maiani et al. [71] predicts a 
mass difference at the level of 8±3 MeV/c^. Recent mea- 
surements of the difference between the A (3872) mass in 
B+ X{3872)K+ and 5° ^ A(3872)A'0 decays [li[l6| 
disfavor this model under the hypothesis that one state 
is dominantly produced in _B+ decays and the other one 
in i?° decays. 

In this Letter we report a study of the mass of the 
A (3872) resonance produced in pp collisions. We con- 
sider the conjecture that the structure observed in our 
data is composed of two different states with distinct 
masses; but failing to discern any evidence for this pos- 
sibility we set an upper limit on the mass difference be- 
tween two hypothetical states. In light of this result we 
perform a precision measurement of the A (3872) mass, 
the main result of this paper. 

The data were collected by the CDF II detector at the 
Fermilab Tevatron pp collider between February 2002 and 
August 2007, and correspond to an integrated luminosity 
of 2.4 fb-i. The CDF II detector ^ consists of a mag- 
netic spectrometer surrounded by electromagnetic and 
hadronic calorimeters and muon detectors. The tracking 
system is immersed in a 1.4 T axial magnetic field and is 
composed of a silicon microstrip detector [lit surrounded 
by an open-cell drift chamber (COT) [l^. It extends out 
to a radius of 138 cm with up to 96 position measure- 
ments in the COT, and achieves a transverse momen- 
tum resolution of (j{pt)/pt ~ 0.15 %pT/(GeV/c). We 
detect muons in planes of multiwire drift chambers and 
scintillators in the pseudorapidity range \rj\ < 1.0. 
Events with J/0 — > decays are recorded using a 

dimuon trigger, which requires two oppositely-charged 
COT tracks matched to muon chamber track segments. 
The reconstructed invariant mass of a dimuon pair is re- 
quired to be between 2.7 and 4.0 GeV/c^. 

To reconstruct X{3872) candidates we first build J/ip 
candidates by combining pairs of oppositely charged 
muon candidates with a transverse momentum, px, larger 
than 1.5 GeV/c. The X(3872) candidates are formed 
by combining J/0 candidates in the invariant mass range 
from 2.95 to 3.25 GeV/c^ with pairs of oppositely charged 
tracks, each with pT > 0.4 GeV/c and assigned the pion 
mass. We require that all four tracks have at least 10 
COT and 2 sihcon hits. For the resulting A (3872) candi- 
dates with pt > 3.5 GeV/c, we perform a kinematic fit in 
which the tracks are constrained to originate from a com- 
mon vertex and the dimuon invariant mass is constrained 
to the world average J/ip mass [2l[. Candidates having 
a kinematic fit of good quality are selected in a broad 
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invariant mass range containing, in addition to X(3872) 
candidates, also ip{2S) candidates that decay to the same 
final state. The '0(25') serves as a valuable control sam- 
ple. 

Several discriminating quantities are combined by a 
neural network into a single selection variable. The in- 
dividual quantities are transformed such that linear de- 
pendences on the invariant mass are removed. The most 
important inputs to the neural network are the Q value 
of the decay, defined as Q = mj/^^+^- — m^+^- — rrij^^, 
the transverse momenta of the two pions, the quality of 
the kinematic fit of the X(3872) candidate, and muon 
identification quantities. The offline muon identifica- 
tion is based on the matching of tracks found in the 
tracking system to track segments in the muon system 
and on the energy deposited in the calorimeter by the 
muon-candidates. For the training of the neural net- 
work, a background sample is extracted from data, se- 
lecting events in regions of the J/07r^7r~ mass away from 
the X(3872) and ip{2S) signals, mainly consisting of J/ip 
particles combined with two random tracks. For the sig- 
nal sample we use simulated X(3872) events. In the 
simulation we generate a single X{3872) per event using 
the momentum distribution of the ■0(25'), which is then 
decayed using the evtgen package Each event is 

then passed through a detector simulation based on the 
GEANT3 package [2j] and a trigger simulation, and is re- 
constructed with the same code as for real data. The 
simulation is in good agreement with the data as veri- 
fied with several kinematic quantities. The final selec- 
tion places a requirement on the neural network output 
and the number of candidates per event. Using wrong- 
sign candidates, where the two pion candidates have the 
same charge, we verify that the selection does not create 
an artificial excess in the mass spectrum. The invari- 
ant mass distribution of the selected candidates in the 
X(3872) mass region is shown in Fig. [TJ The sample 
contains about 6000 X{3872) signal events. 

Before we perform a mass measurement, we test 
whether the signal is consistent with a single state or 
we have evidence for more than one state. In the test 
we perform a binned maximum likelihood fit to the mass 
distribution in data, where we describe the combinatorial 
background by a second order polynomial, and the sig- 
nal by a nonrelativistic Breit-Wigner function convolved 
with a resolution function determined from simulated 
events and parametrized by the sum of two Gaussians. 
The core Gaussian, with a width of 3.2 MeV/c^, accounts 
for two thirds of the resolution function; the second Gaus- 
sian has about twice the width. In the fit we fix the width 
of the Breit-Wigner function to F = 1.34 MeV/c^, our av- 
erage of the widths measured in J/07r"'"7r~ decays (ll[l5|. 
The uncertainty on F of 0.64 MeV/c^ is taken into ac- 
count in the hypothesis test described below. As a test 
statistic we introduce a factor t that scales the intrin- 
sic and resolution widths of the signal shape. The value 
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FIG. 1: Invariant mass distribution of the X{3872) candi- 
dates. The points show the data distribution, the full line is 
the projection of the unbinned maximum-likelihood fit, and 
the dashed line corresponds to the background part of the 
fit. The inset shows an enlargement of the region around the 
X(3872) peak. Residuals of the data with respect to the fit 
are displayed below the mass plot. 



of t determined by the fit to the data is then compared 
to the distribution of t from an ensemble of simulated 
experiments that assume a single state. Based on this 
comparison the consistency of the data with the single- 
state hypothesis is evaluated. The pseudoexpcriments 
are generated using the same fit model as in data. As 
several quantities are known only with limited precision, 
we vary those in the sample generation according to their 
uncertainties. The varied parameters include background 
shape parameters, the number of signal and background 
events, the width of the Breit-Wigner function, and the 
overall width of the resolution function. From a compar- 
ison of the ■0(25) signal in the data to that of simulated 
events we observe that the simulation underestimates the 
resolution by about 5%. The samples were generated 
with a resolution corrected for this discrepancy. 

From data we obtain a width scale parameter value of 
t = 1.052. In Fig. [5] we show a comparison of the fitted 
scale parameter to the distribution obtained from simu- 
lated experiments assuming a single state. We conclude 
that the data are fully consistent with a single state. In 
the absence of evidence for two distinct states we set 
an upper limit on the possible mass difference between 
two hypothetical states. As a test statistic we use the 
width scale t, which is compared to expectations from 
samples simulated with different mass splittings. We as- 
sume that both states have the same mass shape and do 
not interfere. We derive upper limits as a function of 
the fraction /i of the lower lying state to the total ob- 
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yields m^(2S) = 3686.03 ±0.02 (stat)MeV/c2. While this 
value is consistent with the world average ip{2S) mass of 



3686.09 ± 0.03 MeV/c^ |2l|, we use the 60 keV/c^ dif- 
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FIG. 2: Distribution of the width scale t for generated ex- 
periments using the single state hypothesis (histogram) . Also 
shown is the measured value from data (vertical line). 
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FIG. 3: The upper limit on the mass difference Am between 
two states as a function of the fraction /i of the yield of the 
lower mass state. 



served signal. The resulting 90% and 95% C.L. upper 
limits arc shown in Fig. [3l For an equal mixture of the 
two contributing states, the limits are Am < 3.2 McV/c^ 
and Am < 3.6 MeV/c^ at 90% and 95% confidence lev- 
els, respectivel y. This result is complementary to other 
measurements |l5l . ITij in that it does not rely on assump- 
tions about the production of the two hypothetical states 
in versus decays, but depends on /i. 

Lacking any indication of dual states we proceed to 
extract the mass of the X(3872) by performing an un- 
binned maximum likelihood fit using the same fit model 
as used in the previous two-state test. In this fit we fix the 
intrinsic width to F = 1.34 MeV/c^ and the resolution 
parameters to their expected values. Free parameters in 
the fit are the mass of the X{3872), the fraction of signal 
events in the sample, a resolution scale factor, and two 
parameters determining the background shape. 

To check the absolute mass scale we use the nearby 
ip{2S) signal in the same J/i/jTr+Tr" invariant mass spec- 
trum. We use the identical fit model as for the X(3872), 
with the exception that the signal shape parameters are 
adjusted to the world average value of F = 0.337 MeV/ 
[2lj for the intrinsic width, and that resolution parame- 
ters are determined from simulated ip{2S) events. The fit 



ference between our measurement and the world average 
value as an estimate of a possible uncertainty due to un- 
certainties both on our measurement and on the world 
average value. 

Since a possible miscalibration of the momentum scale 
would show up as a dependence of the measured mass on 
momentum, we measure the ip{2S) mass as a function of 
several kinematic variables. We find that any tested de- 
pendence has an effect below 0.1 MeV/c^, which is taken 
as an additional measure of the systematic uncertainty. 
This uncertainty is summed in quadrature with the sys- 
tematic uncertainty on the absolute mass scale derived 
above. To translate the estimation of the mass-scale un- 
certainty from the ip{2S) to the X{3872) we scale the sum 
by a factor of 1.6 that is modeled by a linear dependence 
on the mass with respect to the J/07r+7r~ threshold. This 
yields a total systematic uncertainty of 0.19 MeV/c^ at- 
tributed to the momentum scale. 

To estimate the effect due to the uncertainties in the 
fit model, we refit the data using alternative models. 
These include the use of a linear function instead of a 
second-order polynomial for the background description, 
a single Gaussian function instead of a non-relativistic 
Breit-Wigner function convolved with double Gaussian 
resolution function for the signal description, and fixing 
the natural width F to zero or to twice the nominal value. 
We also perform a fit in a mass window reduced by 40%. 
All of these modifications have a negligible effect on the 
fitted mass, below 20 keV/c^, and therefore we do not as- 
sign any systematic uncertainty to the measurement due 
to the fit model. Because the observed decays to D^D*^ 
may stem from a different particle we assume that the 
mass line shape is not distorted by them. If this were the 
case, as discussed in Ref. [i^l, it would be expected to 
increase the measured mass by about 150 keV/c^. 

The final mass measurement for the X(3872) is 
3871.61 ± 0.16 (stat) ± 0.19 (syst) MeV/c^. The mea- 
sured value is in good agreement with the world average 
[2H and the more precise average of measurements in the 
J/i/'TT+TT" channel including the preliminary Belle mea- 
surement 16[ . It is the most precise single measurement 



to date and improves the precision of the latter average 
by about a factor of 1.5. 

Our measurement is below the D^D*^ mass threshold 
of 3871.80 ± 0.35 MeV/c^ ^ by 0.19 ± 0.43 MeV/c^. 
This implies that the interpretation of the A" (3872) as 
£)0£)*o niolecule is still possible, although the current 
precision does not preclude an ^(3872) mass above the 
]jOjj*o jnasg threshold. A future increase in precision of 
this comparison will therefore require improvements in 
the precision of the D'-* and D*'^ masses. Concerning the 
four-quark hypothesis, our mass splitting upper limits for 
two hypothetical states with relative fractions between 
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0.2 and 0.8 exclude the rane; e of 8 ± 3 MeV/c^ predicted 
in Rcf. 0. 

In summary, we present a new measurement of the 
X{3872) mass using its decay to J/^tt+tt". Our mea- 
sured value of 3871.61 ± 0.16 (stat) ± 0.19 (syst) MeV/c^ 
supersedes that of Ref. 0], and is more than two times 
more precise than the best single measurement so far. 
In addition, we derive upper limits on the mass differ- 
ence for the hypothesis of two X(3872) states, which are 
predicted by some four-quark scenarios, as a function of 
their relative contribution to the observed signal. For 
an equal mixture of the two possible states, the limit is 
Am < 3.6 MeV/c^ at 95% confidence level. 
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